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We extend our analysis of two electrons on a sphere [Phys. Rev. A 79, 062517 (2009); Phys. Rev. 
Lett. 103, 123008 (2009)] to electrons on concentric spheres with different radii. The strengths 
and weaknesses of several electronic structure models are analyzed, ranging from the mean-field 
approximation (restricted and unrestricted Hartree-Fock solutions) to configuration interaction ex- 
pansion, leading to near-exact wave functions and energies. The M0ller-Plesset energy corrections 
(up to third-order) and the asymptotic expansion for the large-spheres regime are also considered. 
We also study the position intracules derived from approximate and exact wave functions. We find 
evidence for the existence of a long-range Coulomb hole in the large-spheres regime, and infer that 
unrestricted Hartree-Fock theory over-localizes the electrons. 
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I. INTRODUCTION 

In recent work, we have reported near-exact [T] and 
exact [2] solutions of the singlet ground state of two elec- 
trons, interacting via a Coulomb potential, but trapped 
on the surface of a sphere. This model was first used 
by Berry and co-workers in the 1980's to provide insight 
into angular correlation in two-electron systems [3H5]. It 
has proven useful for understanding the electronic po- 
larity of nanoclusters and for explaining the giant po- 
larizability of Nai4Fi3 and spontaneous dipole formation 
on niobium clusters [S]. Within the adiabatic connec- 
tion in density functional theory (DFT) [THI] , Seidl care- 
fully studied this system [Tni [TT] in order to test the ISI 
(interaction-strength interpolation) model [H], deriving 
values of the energy by numerical integration. Further- 
more, it has been shown that this kind of spherical con- 
straint applied to the Moshinsky atom leads to a 
solvable Schrodinger equation [T?]. 

Berry and collaborators also considered an extension in 
which each particle is confined to a different, concentric 
sphere [TS] and used this model to simulate the rovibra- 
tional spectra of the water molecule in both the ground 
[l6] and excited states [17]. More recently, the model 
has been applied to quantum-mechanical calculations of 
large-amplitude light-atom dynamics in polyatomic hy- 
drides [iHiiin]. 

It seems timely therefore to generalize our earlier stud- 
ies [2 12] to the case of two electrons located on the sur- 
face of two concentric spheres of different radii. To be 
consistent with our previous work jU 0, we will focus 
on the singlet ground state, which allows us to confine 
our attention to the symmetric spatial part of the wave 
function and ignore the spin coordinates. However, when 



the two radii are not equal, the spin coordinates are ir- 
relevant, and one can easily generalize the present results 
to the triplet state by antisymmetrizing the spatial wave 
function. 

Symmetric and asymmetric Hartree-Fock (HF) solu- 
tions are discussed in Section [III| and the strengths and 
weaknesses of M0ller-Plesset (MP) perturbation theory 
[20] in Section IV We consider asymptotic solutions for 



the large-spheres regime in Section |V] and, in Section jVH 
we study the convergence behavior of the variational con- 
figuration interaction (CI) scheme. Finally, by investi- 
gating the shape of the position intracule and the cor- 
responding Coulomb hole (Section VII I, we report the 
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existence of a secondary Coulomb hole, shedding light 
on long-range correlation effects in two-electron systems. 
Atomic units are used throughout. 



II. HAMILTONIAN 

Our model consists of two concentric spheres of radii 
^1 < -^2) each bearing one electron. The position of the 
i-th electron is defined by the spherical angles (6*^,0^), 
the interelectronic angle 9 by 

cos 6* = cos 01 cos 6*2 + sin^i sin 02 cos(0i — ^2), (1) 

and the interelectronic distance by u = |ri — where 

R2- Ri<u<R2 + Ri. (2) 

The Hamiltonian of the system H is simply 

H^f + u-\ (3) 

where f = fi+f2 = -(Vf V^)/2 is the kinetic energy 
operator and is the Coulomb operator. It is some- 
times convenient to recast H in term of the interelec- 
tronic angle 9. Introducing the dimensionless parameter 
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FIG. 2. R2 X (dotted), i?2 x E^^^ (dashed) and R2 x 

^cxact (gojjfj) as a function of R2 for i?i = 2, 5, 10, 15. 



0<A = i?i/i?2<l and using ([T]), one finds 



2i?2 sin 61 <i6i V dB 



R2VI + A2 - 2Acos( 



(4) 

which shows the different scahng behavior of the kinetic 
and electrostatic terms. 



III. HARTREE-FOCK APPROXIMATION 

A. Symmetric solution 

For i?i = i?2 = R, the restricted Hartree-Fock (HF) 
wave function and energy take [1] the simple forms 



HF 



1 



47ri?2' 



E 



HF 



1 

R' 



(5) 



For Ri < i?2, the electrons occupy different orbitals 
and an unrestricted HF treatment is required. However, 



the high symmetry of the system implies that there is a 
solution in which each orbital is constant over its sphere, 
and the resulting wave function and energy are 



SHF 



1 



47ri?ii?2 



■SHF 



1 



(6) 



We call this the "symmetric Hartree-Fock" (SHF) solu- 
tion, because the orbitals are spherically symmetric and 
we note that the SHF energy depends only on the radius 
of the larger sphere. 



B. Asymmetric solution 

For certain values of Ri and R2, a second, lower-energy 
HF solution arises [2T1123| . in which the two electrons 
tend to localize on opposite sides of the spheres. We call 
this the "asymmetric Hartree-Fock" (AHF) solution for 
the orbitals possess cylindrical, not spherical, symmetry. 

To obtain the AHF wave function 

$^"^(ei,^2) = *l(ei)*2(e2), (7) 

we expand the orbitals as 

00 CX3 

*l(^l) -I]c,*,(^?l), *2(02) =^d£*K^2), (8) 

e=o 1=0 
in the basis of zonal harmonics |Mj 

- Ye{0,)/R, = Y,o{9,,^,)/R,. (9) 
The Fock matrix elements for the two orbitals are 



^1 — 



2Rl 



^ ^3,^4=0 



^ d,,d,,G[% (10) 
,,£4=0 

00 

J2 W,G',% (11) 



where Sg-^^g^ is the Kronecker symbol and 

4 E^TTT^'^^i^^^^^^^^^) (12) 

2 e 

are the two-electron integrals expressed in terms of the 
Wigner 3j-symbols [IS] 

./ (2^1 + 1) (2^ + 1) (2^3 + 1) fei £2 

(13) 

The summation in (12) runs from max(|£i — €2!, — ^4|) 
to min(fi +£2,^3 + because of selection rules [25]. 
The AHF energy is 



00 

^AHF ^ i y 
2 ^ 



^(^+1) 2 ^(^+1) 



(14) 
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FIG. 3. MP2 (dashed), MPS (dotted) and exact (solid) cor- 
relation energies (relative to SHF) as a function of R2 for 
Ri =0.1,0.5,1. 



For all of the radii considered, truncating the expansions 
in (|8| at ^ = 15 yields E^^^ with an accuracy of 10~^^. 
The asymptotic limits of the HF energies satisfy 



lim R2E' 



SHF 



lim {Ri + R2) E^^^ = I, 

Ri .J?2— ^00 



(15) 
(16) 



and the limiting AHF energy corresponds to the Coulomb 
interaction between two electrons that are fully localized 
on opposite side of their respective spheres. Such sys- 
tems are known as Wigner molecules [26j and have been 
observed in a variety of similar systems [57H5D] . 

The creation of localized orbitals leads to decreased 
Coulombic repulsion but increased kinetic energy, and 
an asymmetric solution therefore exists only when the 
former outweighs the latter. By considering an orbital 
basis consisting of only Yq and Yi, it can be shown that 
this occurs only when i?i > i?"" and R2 < 
where i?'^''* — 3/2. Figure [l] illustrates this graphically. 

Figure [2] shows the SHF, AHF and exact energies as 
functions of R2 for several values of Ri. The difference 
^AHF_^cxact (jecpgaggg as Ri = i?2 increases, indicating 
that the AHF energy is asymptotically correct. 



IV. EXPANSION FOR SMALL SPHERES 
A. First-order wave function 

In M0ller-Plesset (MP) perturbation theory gD], the 
total Hamiltonian is partitioned into a zeroth-order 
Hamiltonian Hq = T and a perturbative correction 
V — . The unperturbed orbitals are spherical har- 
monics on each sphere and therefore, from Section p^II A[ 
we have = and E'^'^^ = R^^ ■ 

The -^-th excited eigenfunction and eigenvalue of Hq 



with S symmetry are (TJ [551 IS] 



*K^) = 7^^P?^Kcos0), 



Ei = t{l + 1) 



47ri?ii?2 
1 



1 



2i?2 2Rl 



(17) 
(18) 



In intermediate normalization, the first-order correction 
to the wave function is 



($0 \u ^1 $/) 



Eq — Eg 
1 

27ri?,i 1 + A2 
where it can be shown that 



Z{cos0,X), 



(19) 



1) 



Peix) 



= 1 + log 2 - log(l -\x + \J\ - 2\x + A2) 
log(l - x) - log(A - X + Vl - 2\x A2) 



A 

This yields the normalized first-order wave function 



(20) 



(21) 



1 



2\Rx 



1=1 



£2(^+1)2(2^-1-1) 



B. Second-order energy 



Using (19 1, one finds that the second-order energy 



2(1 - A)2 ln(l - A) + 2(1 + A)2 ln(l + A) - GA^ 



A^ 



H-A2 

2Fi(l,2,5,A) + 2i^i(l,2,5,-A) 



1 + A2 
A^ 4A6 



6 



0(A8 



(22) 



(where 2^'i is the Gauss hypergeometric function [53]) 
depends only on the ratio of the radii. 

When the radii are equal, E^"^^ takes the value 



lim £;(2) =41n2-3, 

A-i-l 



(23) 



which has been discussed by Seidl [TOl [H] and us [H . 
When the radii are very different (i.e. A « 0), the HF 
treatment is accurate and the second-order energy is 

^ CilR\, (24) 

where C4 = — i?f/3. Although ( [24| can be identified as 
the dispersion energy, it does not exhibit the usual R~^ 
behavior. Analogous results have also been reported for 
other systems [33]. 
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C. Third-order energy 



Using (19 1, one finds the third-order energy 



£;(3) = ($(1) 
4Ai?2 



$(1) 



(1 + A2)2 



24A3 + 2A(1- A2)Li2(A2) 
10A(1 - \f ln(l - A) - (1 - A)^ ln2(l - A) 
10A(1 + \f ln(l + A) + (1 + A)-'' ln2(l + A) 



^R2 



2A8 12A^° 
~9 35~ 



0(Ai2) 



(25) 



where Li2 is the dilogarithm function |34j . 

When the radii are equal, E''-^-' takes the value 

lim E^^^ ^ 8i?2(3 - 5 In 2 + In^ 2) 



(26) 



that we have given previously [T] . When the radii are very 
different and i?2 is not too large, E'^^'' is tiny and E'^'^^ is 
a good approximation to the total correlation energy. 
The MP2 and MPS correlation energies, defined by 



E 



MPri 



Y^E^-\ 

m=2 



(27) 



are shown in Fig. [3| For Ri = 0.1, the MP2 and MP3 
energies are accurate for all i?2- For larger i?i, the dis- 
crepancy between the MP and exact energies is noticeable 
for small i?2, but remains small for large R2. 

The MPS energy is usually better than the MP2 en- 
ergy. However, as we have shown previously [1], the MP 
expansion appears to diverge when the radii are similar 
and not small piUfT^ . 



V. EXPANSION FOR LARGE SPHERES 
A. Harmonic approximation 

In the large-spheres (LS) regime, the electrons reduce 
their Coulomb repulsion by localizing on opposite sides 
of their spheres, oscillating around their equilibrium po- 
sitions with angular frequency uj (zero-point oscillations) . 
The same phenomenon has been observed by Seidl and 
collaborators [lOl [El |35H39] . 

In this case, the supplementary angle ^ — n—d becomes 
the natural coordinate of the system. Using the Taylor 
expansions, cotf = k + 0{£,) and 



1 



1 



A 



Vl + A2 + 2Acose (1 + ^) 2(1 + A) 



o(e'), 

(28) 



a> 
c 



re 

o 
o 



0.00 



-0.02 



-0.04 



-0.06 



-0.08 




FIG. 4. E^^° (daslied) and E^^'^ (dotted) correlation energies 
(relative to SHF) as a function of R2 for Ri — 3, 5, 10. The 
exact correlation energy (solid) is also shown. 



the Hamiltonian Q becomes 



JJU, 



2Ri 



A 



(1 + X)Ri 



1 



2{1 + Xy 



The lowest eigenfunction of (|29|) is 
$"(0 oc exp 



yRi 



(1 + A)3(1 + A2) 

and the associated eigenvalue is 

E'^ = 



^LSO -^e — e 



1 



UJ 

2 ■ 



(29) 



(30) 



(31) 



The first term of (31) represents the classical interaction 



of two electrons separated by a distance R1+R2, and the 
second one is the energy associated with the zero-point 
oscillations of angular frequency 



2V(1 + A2)/A 

(i?l+i?2)3/2 • 



(32) 



B. First anharmonic correction 

The first anharmonic correction 

(1 + A^) d _ A^(A2-7A + 1) 4 
QRl ^di 2A{l + \fRi ^ 



W 



(33) 



arises from the next two terms of the Taylor expansion of 
cot^ and the Coulomb operator (28). Defining E^^"^ = 
^Lso _|_ £(1)^ anharmonic correction energy is 



f (1) = A-k'^rIb? 



00 



2, $"(^)p^$-(^)^d^ 


I \2^ 



(1-A + A^)(1 + A^) 
4(i?i+i?2)2 ■ 



(34) 
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The LSO and LSI correlation energies are shown in Fig. 
|4]with respect to i?2 for three values of For the large 
values of both curves agree very well with the exact 
correlation energies, while for the smaller values of the 
radius of the first sphere, LSI systematically improves 
the results compared to LSO. 



VI. CONFIGURATION INTERACTION 

To obtain an accurate wave function, we expand it in 
the Legendre basis 



grows linearly over the domain of allowed u values. The 
AHF intracule is more complicated but is given by 

oo oo 



1=0 



2=0 



with di = V2i + 1 Q and x = (Rl + Rl- u^)/{2RiR2) 
We define the Coulomb hole [HI 



d.V{u)^V{u)-V'^''{u) 



(41) 



(35) 



1=0 



as the difference between the intracule from a correlated 
wave function and that from the lowest HF wave func- 
tion. To quantify the correlation effects, it is useful to 
identify the minimum ii^'^ , the root iP'' and the strength 



where is the CI amplitude of the excited configuration 
The elements of the CI matrix are given by 



\IS.V{u)\du 



(42) 



H 



1) 



h + l) 



2Rl 



2Rl 



1 

R2 



E 

(36) 



of /tJieilho: 



,sr) Coulomb hole. In certain cases, a 
peccMdary long-range (Ir) Coulomb hole appears [32] • Its 
strength is given by 



where (^i £2 f) is given by ( 13 ). 

In our earlier work on the Ri — R2 case [1], we found 
that the CI expansion converges slowly with respect to L 
because of the interelectronic cusp that arises wherever 
the electrons meet [40! . We also showed that this problem 
can be overcome by expanding the wave function as a 
polynomial in u. 

Here, however, we find (Table ^ that the CI expan- 
sion converges rapidly, provided that R2 is significantly 
greater than This is to be expected, because the 
fact that the electrons are confined to different spheres 
means that they can never meet and that the exact wave 
function is therefore cuspless. 



VII. INTRACULES AND HOLES 

To study the relative positions of the electrons in space, 
we have computed the position intracule 



(43) 



V{u)^{<^\5{\ri-r2\-u)\^) 



(37) 



the probability density for the interelectronic separation 
u, from several of the wave functions $ above. Be- 
cause the SHF, MPl, LR and CI wave functions depend 
only on u (or, equivalently, on the interelectronic angle), 
their position intracules are given by the simple Jacobian- 
weighted density 

7'(u) = 87r2i?ii?2ii |$(u)|^ (38) 

For i?i = i?2, the MP2 intracule is also available [T]. 
The SHF intracule 



,SHF 



(«)= 



R2-Ri<u<Ri 
otherwise. 



i?2 



(39) 



5''-= / \Ariu)\du, 
where v}'^ is the long-range root. 

A. Weakly correlated regime 



In the weak interaction limit, the only HF solution 
is the symmetric one (Section III A I and correlation ef- 
fects are well-described by the MP approximation (Sec- 
tion |lV|. 

Figures [5ja) , [Sjjc) and [5]Je) show the Coulomb holes 
derived from the MPl, MP2 and CI wave functions for 
three small values of R = Ri = i?2. For such radii, the 
MP-based and exact position intracules are very similar, 
and become identical as i? — 0. The holes are negative 
for small u and positive for larger u, implying that corre- 
lation decreases the likelihood of finding the two electrons 
close together and increases the probability of their being 
far apart [H]. To illustrate the spatial distributions of 
the electrons, we have plotted the MPl Coulomb holes 
on the surface of a sphere (Fig. |6| for the three same 
values of the R = Ri = i?2 ■ 

The evolution of AP{u) with respect to the increase of 
R2 is shown in Figs, [sjb) , [gd) andjsff). As R2 mcreases, 
the difference between the MPl and exact holes decreases 
and they match perfectly as i?2 — ^ 00. 

Table |TT] shows that the first-order correction reduces 
the probability of small u values too much, and that the 
second-order correction partly corrects this, at least for 
small of values of i?. As a consequence, the strength 
of the MPl hole is always larger than the true one, but 
exhibits the right asymptotic behavior for small R. 



6 



TABLE I. Convergence of correlation energy with respect to the number L of terms in the CI wave function. 



L 


Ri = 


1 i?2 = 


1 


Ri 


= 1 R2 = 


1.1 


Rl^l i?2 = 


1.5 


-Ri = 


1 i?2 = 


2 


1 


-0.131 


665 


623 


696 


-0 


102 


135 


552 


400 


-0.041 049 


324 


810 


-0.015 


832 


811 


848 


2 


-0.141 


241 


198 


782 


-0 


108 


514 


851 


797 


-0.042 647 


613 


578 


-0.016 


238 


558 


022 


3 


-0.144 


065 


402 


167 


-0 


110 


102 


786 


034 


-0.042 870 


984 


637 


-0.016 


271 


251 


778 


4 


-0.145 


273 


783 


726 


-0 


110 


674 


429 


216 


-0.042 915 


528 


628 


-0.016 


274 


960 


819 


5 


-0.145 


900 


461 


200 


-0 


110 


923 


205 


726 


-0.042 926 


179 


125 


-0.016 


275 


462 


730 


10 


-0.146 


847 


645 


782 


-0 


111 


180 


386 


287 


-0.042 930 


208 


158 


-0.016 


275 


553 


424 


15 


-0.147 


047 


095 


403 


-0 


111 


201 


118 


465 


-0.042 930 


221 


857 


-0.016 


275 


553 


441 


20 


-0.147 


120 


296 


106 


-0 


111 


204 


056 


487 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


25 


-0.147 


155 


035 


738 


-0 


111 


204 


595 


556 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


30 


-0.147 


174 


201 


368 


-0 


111 


204 


710 


528 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


35 


-0.147 


185 


880 


267 


-0 


111 


204 


737 


604 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


40 


-0.147 


193 


518 


573 


-0 


111 


204 


744 


445 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


45 


-0.147 


198 


785 


870 


-0 


111 


204 


746 


267 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 


50 


-0.147 


202 


570 


742 


-0 


111 


204 


746 


773 


-0.042 930 


221 


942 


-0.016 


275 


553 


441 



Exact -0.147 218 934 944 -0.111 204 746 979 -0.042 930 221 942 -0.016 275 553 441 



TABLE IL 


Minimum (u) 


root (m) and strength (S) of the 


TABLE 


III. Minimum (u), maximum (u) 


root (u) and 


MPl, MP2 and exact Coulomb holes for _Ri = R2 


= R. 


strength (S) of the exact Coulomb hole for various R. 


R 




Minimum if^ 




R 


Short-range Coulomb hole 




MPl 


MP2 


Exact 




Minimum u'"' Root m'"' 


Strength 5" 


0.1 


0.062 


0.063 


0.063 


1.6 


1.115 2.104 


0.0683 


0.2 


0.127 


0.130 


0.129 


1.7 


0.945 1.836 


0.0148 


0.5 


0.337 


0.353 


0.347 


1.8 


0.787 1.511 


0.0126 


1.0 


0.746 


0.810 


0.757 


2 


0.572 1.081 


0.00295 


1.5 


1.228 


1.413 


1.211 


5 










Root u"' 




10 








MPl 


MP2 


Exact 


100 






0.1 


0.130 


0.131 


0.131 




Long-range Coulomb hole 


0.2 


0.262 


0.264 


0.264 




Maximum v}' Root u^' 


Strength S^' 


0.5 


0.667 


0.678 


0.675 


1.6 






1.0 


1.371 


1.412 


1.386 


1.7 


2.719 3.221 


0.00334 


1.5 


2.109 


2.216 


2.121 


1.8 


2.606 3.151 


0.0267 






Strength 




2 


2.737 3.382 


0.0653 




MPl 


MP2 


Exact 


5 


7.398 8.697 


0.161 


0.1 


0.0245 


0.0235 


0.0235 


10 


15.90 17.95 


0.158 


0.2 


0.048 


0.045 


0.045 


100 


184.62 192.31 


0.134 


0.5 


0.116 


0.096 


0.099 








1.0 


0.211 


0.145 


0.164 








1.5 


0.476 


0.235 


0.210 


For R 


> 1.7, a secondary Coulomb hole appears in 



B. Strongly correlated regime 

In the strong interaction limit, the Coulomb repulsion 
dominates the kinetic energy, an AHF solution exists 
(Section IIIB) and the electrons oscillate around their 



equilibrium positions (Section (|V|). 



creases the probability of finding electrons at large sep- 
arations. This implies that the AHF wave function over- 
localizes the electrons on opposite sides of their spheres, 
and that correlation then delocalizes them slightly. Such 
secondary Coulomb holes are not peculiar to our system; 
they have also recently been observed in the He atom 
and the H2 molecule [43| . 

For R > 2, the primary Coulomb hole disappears com- 
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FIG. 5. On the left: MPl (dashed), MP2 (dotted) and exact (sohd) holes for various iJi = 7?2. On the right: MPl (dashed) 
and exact (solid) holes for various R\ < i?2. 



pletely, leaving only the secondary one (Table III and 
Fig. [7|. The secondary:primary strength ratio is larger 
than in the He atom and the equilibrium H2 molecule 
(1-2%) and resembles that in the H2 molecule at a bond 
length of 3 a.u. [43] . 

Figure [8] shows the evolution of the exact hole for i?i = 
10 and i?2 ranging from 10 to 100. The secondary hole 



vanishes when R2 exceeds Ri/R'^"*' and the AHF solution 
collapses to the SHF one. 

To compare the holes based on the LS wave function 
(Eq. ( 30 )) and the exact one, we have plotted the dif- 
ference between the exact and the LS holes (A AT' (it) in 
Fig. [9]). For R = Ri = R2, the agreement between 
the two holes is fairly good for large R (Fig. |9][a)). 
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R = 0.5 R = 1 R= 1.5 




FIG. 6. MPl holes for various R = Ri = R2 plotted on 
the surface of a sphere. The holes with respect to 9 are also 
represented. 
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FIG. 7. Holes based on the exact wave function for various 
{R = 5, 10, 100). The inset graph shows holes for radii just 
above the critical value Tl'^"* = 3/2. 
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FIG. 8. Holes based on the exact wave function for _Ri = 10 
and various R2 (from 10 to 90). The inset graph represents the 
transition between the AHF and SHF solution {R2 = 200/3). 
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FIG. 9. At the top: difference between the exact and the LS 
holes (AAP(it)) for various _Ri = R2: 5 (solid), 10 (dashed) 
and 50 (dotted). At the bottom: AAV{u) for Ri = 10 and 
various R2 (10, 30 and 50). 



For the smaller values of the radius, it shows that the 
electronic zero-point oscillations tend to over-localize the 
electrons compared to the exact treatment. However, 
the secondary Coulomb hole is less pronounced but still 
present in the LS approximation. Moreover, one can see 
that the LS treatment slightly increases the likelihood of 
finding the two electrons close together. 

Figure ^h) reports the modification of AAV{u) for 
a fixed value of the first sphere radius = 10) and 
various R2 (10, 30 and 50). When R2 is increasing, the 
first minimum disappears, and the main effect of the LS 
approximation is thus to over-localize the electrons on 
opposite side of the spheres. 

The 2D spatial distribution of the electrons is depicted 
in Fig. [TOj where we have represented the LS holes for 
various R = Ri ~ R2 (5, 10 and 50) on the surface of a 
sphere. 

VIII. CONCLUSION 

We have performed a comprehensive study of the sin- 
glet ground state of two electrons on the surface of 
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R = 5 R = 10 R= 50 




FIG. 10. LS holes for various R = Ri — R2 plotted on the 
surface of a sphere. The holes with respect to 6 are also 
represented. 



totic solutions of this problem, we have taken into ac- 
count, in the harmonic and anharmonic approximations, 
the zero-point oscillations of the electrons around their 
equilibrium position. For any values of R2 > Ri, the 
near-exact wave function and energy can easily be ob- 
tained by a CI expansion based on Legendre polynomials 
because there is no cusp in the wave function. 

A study of the position intracules and Coulomb holes 
reveals that, as in the helium atom and hydrogen 
molecule, there is a secondary Coulomb hole in the large- 
sphere regime. Indeed, as R increases, the primary hole 
disappears and only the secondary one remains. This 
reflects an overlocalization of the electrons in the asym- 
metric Hartree-Fock solution. 

Our results should be useful for the future develop- 
ment of accurate correlation functionals within density- 
functional theory [TTl [38l [39l |44] and intracule functional 
theory [45H50] and, also, for understanding secondary 
Coulomb holes in more complex systems [351 SS] ■ 



spheres of radius i?i and R2 ■ Symmetric and asymmetric 
HF solutions show that the symmetry-breaking process 
occurs only when i?i > i?™* = 3/2 and R2 < Rj/R"'^. 
MP2 and MPS energy corrections reveal that MP theory 
is appropriate when both radii are small (as previously 
known) and, also, when R2 3> Ri- To derive asymp- 
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